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r~| ' Abstract - The qualitative analysis of the initial value problem V related to a non linear 

C^ ' third order parabolic equation typical of diffusive models is discussed. Some basic properties 

of the the fundamental solution of a related linear operator are determined and are applied 
to an equivalent integro differential formulation of the problem. By the fixed point theorem, 
existence and uniqueness results are obtained. 
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^ : 1 Introduction 



f\ • Several mathematical models in applied sciences - as Biochemistry, Epidemiology, 

c^ , Population dynamics - deal with higher order parabolic equations typical of diffusive 

phenomena. The results of experimental analysis show that the diffusion coefficients 
are generally small and this requires too much time in the transmission of signals. As 
consequence it necessary to correct the diffusive models by means of reaction-diffusion 
systems. So, one obtains non linear hyperbolic equations perturbed by viscous terms 
of higher order. 

A typical example is given by equations as 

(1.1) Lu= {edt + c^)u^x- {dt + a)ut= F{x,t,u,Ux), 



which are object of an extensive hterature [1]- [5]. 

When F is hnear, numerous contributions are already known [6] - [15]. In partic- 
ular the fundamental solution K of the operator L has been obtained in [16] where 
some basic properties has been discussed too. 

In this paper, the non linear case of (1.1) will be considered and existence and 
uniqueness results for the initial value problem V in all of the space will be obtained. 

For this, the equivalence between the problem V and an integro - differential 
equation with kernel K is deduced and the fixed point theorem is applied by means 
of appropriate estimates for K. 

This analysis can be also applied to obtain the dependence continuously upon the 
data, stability properties and boundary - layer estimates. More, the study of other 
boundary value problems (semi infinite or finite media) with Dirichlet or Neumann 
conditions could be achieved too. 



2 Fundamental solution and its properties 



If D = {{x,t) : X E ^ , < t <T }, let r = \x\/e, b = c^/e. As one can verify, the 
fundamental solution of the L— operator is simbolically defined by 



(2.1) K{r,s) 
If one puts 

(2.2) G(r,s) 
obviously one has: 
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(2.3) k{r,s)= -^G{r,s) => K {r,t) = f G{r,'. , 

VS -^0 y/7r(t - r) 

where G{r,t) is the inverse tranform of G{r,s) given by (2.2). 



If In{z) denotes the modified Bessel function of the first kind, the exphcit expres- 
sion of G{r,t) is given by : 



(2.4) G{r,t) = —^ f ^ e-^(*"^) i^(r^(b - a){t - v)/v dv, 
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and the following theorem holds: 

Theorem 2.1- For all r > 0, the Laplace integral Ct G{r,t) converges absolutely 
in the half-plane ^es > max{ —a, —b), and one has: 

p - r ^/s{s+a)/(s+b) 

(2.5) CtG{r,t) = G{r,s) 



2^^{s + a){s + b) 



According to the results extablished in [16], for the function K the following prop- 
erties hold : 

Theorem 2.2- The function K defined by (2.3), (2.4) is a G°^{D) solution of 
the equation Lu = 0. When a < b = c^/e, K is never negative in D and more it 
results: 

/oo 1 

K{\x-ilt)di = (-) (1-e-"*) < 1/a, 
-oo Oj 



(2.7) lim drK{r,t) 
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Moreover, one has: 

Corollary 1.2 - By means of the Laplace transforms, it's easy to verify the fol- 
lowing formulae: 

(2.8) / (dt +a)K (x,t) dx = 1, 
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(2.9) / (dt +b)K (x,t) dx = e"'^* + 6(1 - e-'^*)/a 

(2.10) / {dt +a - edl^)K {x,t) dx = 1 - e"^* . 



3 Initial value problem and explicit solution in the 
linear case 

Let consider now the initial value problem V in D : 

Lu= {edt + c^ ) u^x - {dt + a)ut = F{x,t,u,Ua:), {x,t) e D, 

(3.1 

u{x, 0) = fo{x), ut{x, 0) = j\{x), X e 3fJ 

where the initial data /o and /i are arbitrary specified functions. 
If g{x) is a continuous function on 3?, consider the convolution: 

(3.2) Ug= f g{OK{x-^,t)d^ = K*g. 

Then, owing to the properties of K{x, t), it's possible to prove the following theorem: 
Theorem 3.1- // a, (3, 7 are three positive constants such that 

(3.3) \g{x)\ < a e^l'l'^' 0<7< 1, 

then the function Ug is a C°°{D) solution of the equation Lu = 0, such that 



(3.4) Imi Ug{x,t) = 0, lim dtUg{x,t) = g{x) 



uniformly for all x in a compact subset of — oo < a; < oo . I 

Observe now that the convolution K * g is infinitely differentiable in D and, 
according to the results of section 2, one has: 



^S_5)^0t+a-edl)u,=e-''g{x)+ f g{^) {dt + a - edl)K{x - ^,t) d^ 
and (2.10) implies: 

(3.6) hm / g{0 (dt+a-edl) K{x - U) d^ = 0. 

Consequently, one has : 

Theorem 3.2 - // g{x) verifies the hypotheses of Theorem 3.1 , then the function 

(3.7) u*g = {dt + a - e dl^) Ug 
represents a smooth solution of the equation Lu = 0, such that 



(3.8) lim u*g{x,t) = g{x), lim dtu*g{x,t) = 

uniformly for all x in a compact subset of — oo < a; < oo . I 

Remark 3.1 - When F = , the results of Theorems 3.1 - 3.2 allow to assert 
that the function 



(3.9) u{x, t) =Uf^ + {dt +a - edl^) Uf, 

represents a smooth solution of the homogeneous initial - value problem (3.1). I 

Consider now the case that F is not vanishing but it is a linear known function 
F = f{x,t) and let Uf the volume potential 

(3.10) Uf = f K*f dT = f dr f f(^,T)K {x-U-r)d^ . 
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It's possible to prove the following 

Theorem 3.3 - When the source term f {x,t) is continuous and everywhere 
bounded function in the set — oo < ^ < oo, < r < T, then Uf represents a smooth 
solution of the problem Luj = — f{x, t) with vanishing initial conditions. 

Proof - If 1 1 / I It = sup 2) \ f {x,t)\, by (2.6) one has: 

(3.11) I M/ (x, t) I < a"^ t 1 1 / I It ^ lim m/ = . 
Further it results: 

(3.12) \K*f\ < a-Ml-e""^*""M II/IIt 
and so 

(3.13) ]imJ^K{x-^,t-T)f{tT,)d^=0 ^]\mdtUf = 

More, by means of estimates similar to (3.6), one deduces that: 

(3.14) hm / f(^,r){dt+a-edl)K{x-C,t-T)d^= 0, 

and this, together with the properties of K (theorem 2.2), implies Luf = — f{x,t). 
Further, theorems 3.1 - 3.2 imply that Uf{x,0) = 0, dtUf{x,0) = 0. I 

At last, by (3.9) and theorem 3.3, it follows: 

Theorem 3.4 - When F = f {x,t) and the data f, /o, /i verify the hypothe- 
ses specified in theorems 3.1 - 3.3, then an explicit smooth solution of the problem 
(3.1) is given by 



(3.15) u{x,t) = —Uf + Uf-^ + {dt + a — e d': 
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4 Integro differential formulation of the non linear 
problem 



Consider now the non linear case of the problem (3.1), where F = F {x,t,u,Ux) is 
defined on the set 



Q = {{x,t,u,p) : {x,t) E D, —oo<u<oo] —oo<p<oo}, 
with D = {{x,t) : X e ^ , < t < T }. 

From now on we shall assume for F = F[x, t, u{x,t), Ux{ x,t) ) the following 
Hypotheses H: 

- The function F{x,t,u,p) is defined and continuous on Q and it's bounded for all 
u and p. 

- For each k > and for \u\, \p\ < k, the function F is Lipschitz continuous in x 
and t for each compact subset of D. 

- There exists a constant f3F such that: 

\F{X, t,Ui,Pi) - F{x, t, M2, P2) I < /3f { |Mi - M2 I + bl - P2 I } 

holds for all (mj, pi) i = 1,2. 

Let M be a smooth solution of the problem (3.1) and, referring to the function 

K{x — C,, t — t) u{C,, r), consider the operators Ai{x, t, C,, r,) , J\f{x, t, C,, r, ) 
defined by: 

^A (x, t, ^, T,) = e K u^^ — K Ur + u Kr — aK u 
(4.1) 

A/'(x, t, ^, T,) = u{ edr — c^ )K^ — u^{ edr — c^ )K 

Further, let L be the adjoint operator of L: 
(4.2) LK = -e K^^r + c^ K^^ - K^r + aK. 



T 



in order to have: 

(4.3) KLu - uLK = KF. 
As consequence, one has 

(4.4) drM +d^M =KF. 
If Up denotes the non hnear potential volume 

(4.5) UF = tdr f K{^-x,t-T) F(e, r, u{^, r), u^{^, r) ) d^ , 

Jo Jn 

it suffices to integrate (4.4) on D in order to have 

/oo 
[M{x,t,i,t) -A^(x, t, e, 0)]rfe =UF 
"OO 

where the term with d^M doesn't contribute because lim|^|_^oo A/" = . 

Let assume that /o G C^ {^), /i G C^ ($R), and more /o, /g, /q', /i let bounded 
on 3fJ. Then, by means of (3.4) - (3.6) and estimates similar to (3. 13), (3. 14), by (4.6) 
one obtains the following representation 

(4.7) u{x,t) = -Up + f /i(0 K{^-x,t) d^ + 



+ (dt + a- 



ed,,) f fo{OKU-x,t) dC 



Viceversa, suppose that the integral differential equations (4.7) has a solution 
u = u{x, t) such that u and Ux are continuous and bounded for a; G 3f? and < t < T. 
Then the terms depending on the initial data /o, /i are differentiable for t > with 
bounded derivatives for t > . The properties of the fundamental solution K and the 
hypotheses H allow to prove that also the volume potential up and its derivatives are 
Lipschitz continuous. Then the function u characterized by (4.7) verifies the problem 
(3.1) everywhere in D. 



The complete equivalence between the problem V and the integro differential 
equation (4.7) confirms the existence and uniqueness of the solution of (3.1) as soon 
as when this will be demostrate for (4.7). 



5 Existence and uniqueness 

Consider now a time - interval [O,?]], with rj < T. At first we will show that the 
solution of the integral differential equation (4.7) exists and is unique for t G [0,1]] 
and then we will extend this for any finite T. 

For ^ > 0, let 

(5.1) Dr^ = {{x,t) -.xe ^, te[0,r]]} 
and let i3^ the space 

(5.2) B^ = {v{x,t):v, v^e CiD^)} 
equipped with the norm 

(5.3) ll'^llr? = sup |f(a;,t)| + sup \vx{x,t)\ < oo. 

The set is a Banach space. The arguments discussed previously allow to state 
that the mapping 

(5.4) J^v{x, t)= f hiOm - ^, t)di + {dt+a- ed^x) f /o(0^(e - x, m+ 

rt roo 

+ dr K{x-U-r)F{^,T,v{^,T),v^{^,T))d^ 

Jo J —oo 



maps Br^ into i3^. 

Estimate now the continuity of J-'. For all t G [0, rj], owing to (2.6) and hypotheses 
H , it results 



(5.5) |J-'fi(a;, t) — J^t>2(x, t)| < /3f a ^ t\\vi — V2\\ri < /3f a "^ ^ H'^^i — "^^2! 

As for the x - derivative, by means of Kr{z,t) , one can prove that: 

d d 1 

(5.6) \ — J'vi{x,t)-—J'v2ix,t)\ < (3f a~^ ^ = r] ||t;i - t;2||r;- 

'^x Ox \ eih — a) 



As consequence, by (5.3), (5. 5), (5. 6) it follows 

(5.7) \\J^vi{x,t) - J^V2{x,t)\\r^ < (3f (- + . =) V \\vi -V2\\n- 

« Jeib - a) 



If we select rj such that 
(5.8) /9f (-+ , ^^ ) <1, 

then J-" is a contraction of i3^ into B^i and so J-" has a unique fixed point u{x, t) E Bf] ■ 

In order to show the existence and uniqueness of the solution in all [0, T] we 
proceede by induction. 

Let assume that the equation (4.7) admits a unique solution u which is bounded for 
< t < kr] together with u^ ■ 

For krj < t < (k + 1)1], let consider the space 



(5.9) Br^ = {v{x,t) : v{x,t-kr]) E 8^} 

and the mapping 

{5.10yF^v{x,t) = f f,(^)K{^-x,t)d^ + {dt + a-ed,,) f fo{OK{^-x,t)d^+ 

/•krj rOD 

+ dr K{x-^,t-r)F{^,r,u{^,T),u,{^,r))d^ + 

Jo J~OD 
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rt roc 

Jkrj J—oo 

which maps Br, into Br^. Moreover, also the mapping J^i satisfies estimates hke (5.7). 

As consequence, for all t G [0, (/c + 1)1]] the mapping J-*! is a contraction and 
hence admits an unique fixed point u & Br,. 

We remark that , as (5.10) shows, the functions u an u^ are continuous also on 
the junction t = Krj. 

In conclusion, the results in sections 4 and 5 allow to state the following theorem: 

Theorem 5.1 - Let the initial data such that /o G C^ ($R), /i G C^ ($R). Let 
/o,/o, /o', /i bounded on ^. Further, let F verify the hypotheses H . Then, the non 
linear problem (3.1) admits a unique regular solution in D. 
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